Tunneling Density of States, Pair Correlation, and Josephson Current in 
Spin-Incoherent Luttinger Liquid-Superconductor Hybrid Systems 
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We study a hybrid system consisting of a spin-incoherent Luttinger liquid adjoined at one or 
both ends to a superconductor. We find the tunneling density of states diverges at low energies 
and exhibits a universal frequency dependence independent of the strength of the interactions in 
the system. We show that in spite of exponentially decaying pair correlations with distance into 
the spin- incoherent Luttinger liquid, the Josephson current remains robust. Compared to the zero 
temperature Luttinger liquid case there is a factor of 2 reduction in the critical current and a halving 
of the period in the phase difference between the superconductors. Our results open the way for a 
new class of experiments in the spin-incoherent regime of one dimensional systems. 
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The low energy behavior of gapless one dimensional 
systems, which may be realized by electrons in quantum 
wires or nanotubes, is described by Luttinger liquid (LL) 
theory The elementary excitations of the interact- 
ing theory are decoupled bosonic charge and spin modes 
that propagate with different velocities, a phenomenon 
known as spin- charge separation that has already been 
observed experimentally Strong repulsive interac- 

tions tend to suppress the spin velocity while enhanc- 
ing the charge one, thereby accentuating the spin-charge 
separation. For strong enough interactions a window of 
energy opens at finite temperature where the spin sector 
consists of thermally excited (randomized) states while 
the charge sector is essentially at zero temperature. A 
one-dimensional (1-d) system in this regime is known as 
a spin-incoherent Luttinger liquid (SILL) 3]. 

While the theory of the SILL has progressed rapidly [3], 
the challenge of reaching the desired window of energy 
has slowed experiment. To date, the best experimental 
evidence has appeared in momentum resolved tunneling 
on gated quantum wires Unfortunately, the analy- 
sis of the experiments is somewhat involved [3j and it 
has become highly desirable to propose (and carry out) 
new experiments to probe the SILL. Recent experimental 
progress has made it possible to fabricate devices consist- 
ing of nanotubes or quantum wires 0] between two 
superconductors (SC). Through gating to modulate the 
electron density and interaction strength, such devices 
open the possibility of studying the SILL in SILL-SC hy- 
brid structures. 

In this Letter we address the theory of SILL-SC hybrid 
structures. Various aspects of LL-SC structures have 
been discussed in the literature already, including the 
tunneling density of states [1| , pair correlations [§| , and 
Josephson current Compared to the LL case, 

the SILL-SC structures exhibit a number of remarkable 
features. In particular, we find that the tunneling den- 



sity of states of a SILL contacted to a SC diverges at 
low frequencies with a universal form independent of the 
strength of the interactions in the system. By contrast, 
for the isolated LL, SILL, and the LL-SC system the 
energy dependence of the tunneling density of states de- 
pends on the strength of the interactions. We also com- 
pute the decay of the pair correlations into the SILL from 
the SC and, as might be expected from the highly ex- 
cited spin states, they decay exponentially fast with a 
length scale set by the interparticle spacing. However, 
the Josephson current resulting from the coherent prop- 
agation of Cooper pairs through a finite length SILL re- 
mains robust. Moreover, the critical current shows the 
same scaling with length as in the case of a LL but its 
value is reduced by a factor of two. As a function of 
phase between the two superconductors, the period of the 
Josephson current is halved. Both the tunneling density 
of states and the Josephson current should be experimen- 
tally accessible. Observation of the results described here 
would be a smoking gun signature of the SILL in these 
hybrid structures. 
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FIG. 1: (color online) Schematic of the model we study. A 
spin-incoherent Luttinger liquid (SILL) of length L is adiabat- 
ically connected to two superconductors (SC) with phase xi 
and X2, and the same superconducting gap A. The tunneling 
density of states a distance x from the end of the SILL can 
be probed with electron tunneling from a metallic lead such 
as scanning tunneling microscope (STM) tip. A Josephson 
current flows through the SILL when x = Xi ~ Xa / 0. 
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A schematic of our model is shown in Fig.l. We assume 
the SILL is adiabatically connected to the SC so that the 
scattering at the interface is in the Andreev hmit. In 
all our calculations we assume that kgT <^ A where kg 
is Boltzmann's constant, T is the temperature, and A 
is the magnitude of the superconducting gap. Bosoniza- 
tion procedures have been developed for LL-SC hybrid 
systems [lH valid for energies w <C A <C Ep that are 
also small compared to the characteristic spin and charge 
energies. At the lowest energy scales the Hamiltonian is 
given by if = Hp + where 



dx 
2^ 



1 

l9p 
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and has the same form only with p replaced by 
a. Here 9p{x) and <j)p{x) £^re bosonic fields representing 
charge and current density fluctuations, gp measures the 
strength of the interactions {gp = 1 for non-interacting 
systems, gp < 1 for repulsive interactions, and — i for 
SU{2) invariant spin interactions), and Vp is the velocity 
of the charge modes. The mode expansions of the bosonic 
fields for a LL between two superconductors separated by 
a distance L are M 
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"■^ with a a short distance cutoff 



where 7, ^ y 

on the scale of the intcrparticlc spacing, [bf3q,b^p, ^i] = 
Spf3'Sqq' for [3 = p^a, and x is the phase difference of the 
order parameter of the two superconductors (assumed to 
have the same A). The integers J' and M are the topo- 
logical (zero mode) numbers that are related to ex- 
cess charge and spin densities. They obey the constraint 
J' + M = even. Substituting the expansion ^ into ([T]) 
leads to 
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While the form of the energy contribution from the non- 
zero modes (g > 0) is valid only for low energies relative 
to the spin and charge energy, the zero mode contribution 
is valid at all energies, in particular in the spin-incoherent 
regime defined by the condition E'spin ksT <^ £^chargo 



spin/charge 



for fiuctuating quantities like the 
Greens function and pair correlations. In terms of the 
fields ([2]) the bosonized electron annihilation operator is 

where s =t, i- The charge fields 9p — {6^ + 6i)/y/2 and 
the spin fields 9a = {9^—9i)/ \f2 with identical definitions 
for 0p and <i>„. 

Single particle Greens function. We compute 
the single particle Green's function Gs{x,x';t,t') = 
-~{Tr'4!s{x,T)tpl{x' ,t')) for a SILL connected to a sin- 
gle SC. Fiete and Balents have developed a simple 
but powerful method to evaluate such expectation values 
and we employ it here. Let us first consider the trace over 
the spin sector, which we assume consists of highly ther- 
malized random spins. The dominant contribution comes 
from the terms where there are effectively no exchanges 
of particles for all the particles between {x' , t') and (x, r). 
For the single particle Greens function, this implies that 
all spins have same orientation between {x' , r') and (x, r) 
and this occurs with a probability of 2~l''^(^''^'^ '"^ )l where 
N(x^ t; x' , t') is the number of electrons between the two 
points. A factor (— 1)^(^''^;^ '"^ ) arises from the permuta- 
tion of the propagating electron with the other electrons 
in the SILL. The general result is [IS] 

00 

GAx,x';t,t')^ ^ {5{m-N{x,T;x',T')){-ir 

m— — 00 

where the remaining expectation value is taken over the 
charge degrees of freedom at zero temperature. The num- 
ber operator is related to the 9p field. 



N{x, r; x',t') = n{x-x') + 



V2, 
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with n = ^ = 2kp/7T, where kp is the Fermi wavevector 
and a is the mean interparticle spacing. For a system 
with a boundary, N and Gs are not space translationally 
invariant [l^ . 

When |a; — a;'| is small, the Green function can be ex- 
pressed as 
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where 0p and ^p are 6p(a;, t; x', r') = 9p{x, T) — 9p{x', t') 
and <I>p(x, t; a;', r') = <pp{x,T) — 4>p(x' ,t'), respectively. 
For X = x' and r' = 0, the expansions with L 
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for zero mode properties and The resulting Green's function (valid for x, VpT > a) is 
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where we have kept only the dominant m = term [l^ . gives 



Tunneling Density of States. The local tunneling den- 
sity of states, A{x,uj), can be computed by Fourier trans- 
forming the Greens function. For VpT ^ x, the Green's 

1 



function Gs{x, x;t,0) 
implies 



(2a-) 



y/gp ln{vpT/a) 
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Thus, the tunneling density of states of a SILL-SC hybrid 
system diverges at low energies in a universal manner in- 
dependent of the interaction strength gp. The result ^ is 
quite dramatic in light of known results for other related 
systems. If a LL with SU (2) symmetric interactions is at- 

tached to a SC, All-sc(2^, l^) - ' ^ ' '--^(sp-i) 



{2xy 



which also diverges for repulsive interactions {gp < 1) but 
in an interaction dependent way 0] . By studying the fre- 
quency dependence of the tunneling density of states for 
different distances from the boundary (with a scanning 
tunneling microscope, for example) both gp and whether 
or not the 1-d system is in the LL or SILL regime can 
be deduced. By fixing a position x and changing the 
density of the wire with a back gate as in Ref.l^], gp 
can be tuned allowing a test of the interaction depen- 
dence (or lack thereof) of the low frequency behavior of 
A{x, Lu) thus enabling an unambiguous determination of 
SILL physics. 

Pair Correlations. We now turn our attention to the 
pair correlation, F{x) = —{ipij{x)'ipi{x)), the most nat- 
ural measure of the proximity effect. We assume the 
Cooper pairs leak into the SILL at x' = 0, and calcu- 
late their amplitude at distance x > a. The calcula- 
tion in the spin-incoherent case closely follows that of 
the Green's function sl, [Hj. The no particle exchange 
condition here tells us that the dominant contribution 
from the spin sector is an alternating up/down configu- 
ration. The probability of finding an up/down configu- 
ration from the boundary to a distance x is 2~l^(^'°'°'''^l 
where N is given by Fermi statistics gives a fac- 

tor (— which trivially evaluates to unity (because 
the Cooper pair is effectively a boson), and the annihila- 
tion of the two electrons is accomplished by the operator 
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so that the pair correlations decay exponentially fast into 
the SILL. Note that the exponential factor is the same 
as that found for the single-particle Greens function in 
an isolated SILL [l^. This is, in fact, a general result 
in the spin-incoherent regime: all n-body particle non- 
conserving operators decay with the same exponential 
factor because the no exchange condition is a result in- 
dependent of n (for sufficiently large x). It is instructive 
to compare the result 1^ with that obtained [^, [3l for 
a LL when T < v^/a <^ Vp/a using the mode expan- 



sions ©: F{x) ^ e 
spin correlation length £,{T) 



' U(T) 



where the 

Va/iTrkBT). As discussed 
in Ref.[16] the SILL can often be understood from the 
LL result when ^(T) saturates to a number of order the 
mean particle spacing: £,{T) = 0{a). Clearly this is also 
the case for the pair correlations, aside from the correc- 
tion — ^ i^ir) exponent in ([9]) from Gaussian 
density fluctuations (reminiscent of the Greens function 

Josephson Current. It is a remarkable result that a cur- 
rent of Cooper pairs flows between two superconductors 
when the phase difference x between them is non-zero 
(See Fig. [1]). The Josephson current may be obtained 
from the well known thermodynamic relation [lo| 
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d\nZ 
dx 
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where Z is the partition function of the 1-d system and 
— e is the charge of the electron. The LL case has been 
investigated before [H, [l3| • Our approach here is to as- 
sume our system is initially in the LL regime at T = 
and then take the Va-/L <C fc^T <ti Vp/L temperature 
limit [l6| using the mode expansions ^ whose zero mode 
components are valid at all energies, in particular in the 
spin-incoherent regime. At all energies the "J' + M — 
even" constraint relates the charge and spin parts of 
the topological (zero mode) terms [9] of the Hamiltonian 
([3]). Since the topological terms are decoupled from the 
non-zero modes, the partition function factorizes, Z = 

Z{x)Z, where Z{x) = Ej'+m- - - 



e *20p(a:)/%/2^ Combining thcsc after taking the Spin tracc is the partition function of the zero modes with ep 
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. "^r? rr and Z describes the non- 
zero modes. To remove the "J' 



connected to Fermi hquid leads, where the conductance 



M = even" con- 
straint, we let J' = 2j + ji and M — 2m + rrii, with 
j,m = 0,±1,... and sum over the two sectors ji — 
m, = 0,1. Thus, Z{x) = Ej.™e-^-(2j+f)%-..(2m)^ ^ 

g-.,(2j + l+f)^g-e„(2m+l)^^ j^O^-g Z + 2tt) = Z{x)\ thuS 

we restrict our domain to |x| < tt. For low temperatures 
and high temperatures relative to both charge and spin 
velocities, J(x) has been worked out in Rcf.[9]. Here we 
focus on the spin- incoherent case where e^. ^ 1 <C Cp. 
Our exact evaluation of J(x) is plotted in Fig. [5] for 
Ccr/ep = 0.1 and €cr/^p = 0.5. 
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FIG. 2: (color online) Josephson current as a function of phase 
difference x between two superconductors and temperature 
in units of Ep = ■KVpgp/{AL) = €p{kBT). Top: E^/Ep = 0.1. 
Bottom: E„/Ep = 0.5. Here E„ = -Kv^/ig^AL) = e^iksT). 

At the lowest temperatures the zero temperature re- 
sult J(x) = (—77^) f coming from the rm — ji = 
terms is well approximated in both cases. However, at 
slightly higher temperatures contributions appear from 
the rrii = ji = 1 terms in the partition function and cre- 
ate an additional zero in J(x) for 7r/2 < |x| < tt. The 
location of this zero depends on temperature, but for 
temperatures low compared to Ep the zero is most sensi- 
tive to the ratio of spin to charge energy, E^/Ep. Thus, 
the zeros of J(x) may be used to deduce the ratio EcjEp. 
In the limit E^jEp <C 1 this zero moves to ±7r/2. As Fig. 
[2] shows, already iox EajEp — Ql the limit iJ^Z-Bp ^ 1 is 
approached. The condition E^/Ep <C 1 is a prerequisite 
for SILL physics, so J(x) can establish whether the 1-d 
system has large enough spin-charge separation for SILL 
physics. In the extreme limit Ep ^ ksT ^ Ec ^ 0, the 
exact form of the spin-incoherent Josephson current can 
be obtained analytically: 
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f < X < 
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Compared to the LL case at T = 0, both the period and 
the critical current are halved in the SILL regime. The 
latter is reminiscent of Matveev's result for the conduc- 
tance of a quantum wire in the SILL regime adiabatically 



reduces to ^ 
value of ^ 



per mode rather than the T 
Thus 



universal 
the observation of a Josephson 
current that follows (fTTj) is a clear indication of SILL be- 
havior. It is remarkable that in spite of the exponentially 
decaying pair correlations ((9]) the Josephson current re- 
mains robust, only with a factor of two reduction in the 
critical current compared to the T — LL result. Phys- 
ically, this is because the SC phase difference x couples 
only to the charge degrees of freedom which remain co- 
herent in the SILL. Note also that the critical current in 
the SILL regime scales as ~ 1/L as general arguments 
rec^uire [9j]. 

In summary, we have determined the properties ex- 
pected for a SILL adiabatically connected to one or two 
superconductors. The tunneling density of states exhibits 
a universal frequency dependence independent of inter- 
actions in the system and the Josephson current has a 
saw-tooth form with a factor of two reduction in the crit- 
ical current and a halving of the period. If the contacts 
of the SILL are non-ideal (non-adiabatic), the adiabatic 
regime may still be obtained at low energies as impurities 
are irrelevant in the SILL for ^p > 1/2 [16|. In this sense, 
the adiabatic model is even more relevant for the SILL 
than for the LL. 
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